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ABSTRACT 

Anomalous microwave emission is known to exist in the Perseus cloud. One of the 
most promising candidates to explain this excess of emission is electric dipole radia- 
tion from rapidly rotating very small dust grains, commonly referred to as spinning 
dust. Photometric data obtained with the Spitzer Space Telescope have been repro- 
cessed and used in conjunction with the dust emission model dustem to characterise 
the properties of the dust within the cloud. This analysis has allowed us to constrain 
spatial variations in the strength of the interstellar radiation field (xisrf), the mass 
abundances of the PAHs and VSGs relative to the BGs (Ypah and Yvsg), the column 
density of hydrogen (Nh) and the equilibrium dust temperature (Tdust)- The param- 
eter maps of YpAH, YysG and Xisrf are the first of their kind to be produced for the 
Perseus cloud, and we used these maps to investigate the physical conditions in which 
anomalous emission is observed. We find that in regions of anomalous emission the 
strength of the ISRF, and consequently the equilibrium temperature of the dust, is 
enhanced while there is no significant variation in the abundances of the PAHs and 
the VSGs or the column density of hydrogen. We interpret these results as an indi- 
cation that the enhancement in Xisrf might be affecting the properties of the small 
stochastically heated dust grains resulting in an increase in the spinning dust emission 
observed at 33 GHz. This is the first time that such an investigation has been per- 
formed, and we believe that this type of analysis creates a new perspective in the field 
of anomalous emission studies, and represents a powerful new tool for constraining 
spinning dust models. 

Key words: ISM: individual(Perseus cloud, IC 348) - radiation mechanisms: gen- 
eral - ISM: abundances - dust, extinction - infrared: ISM 
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1 INTRODUCTION 

Cosmic Microwave Background (CMB) anisotropy measure- 
ments have provided a great insight into the cosmological 
parameters that define our Universe. Obtaining these mea- 
surements to ever higher sensitivity is complicated by the 
presence of contaminating foregrounds, whose physical un- 
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derstanding is therefore critical. The discovery of a dust- 
correlated emission mechanism in the frequency range 10 - 
100 GHz (Kogut et al. 1996; Leitch et al. 1997; de Oliveira- 
Costa et al. 1997) has reignited the study of Galactic fore- 
grounds as interstellar medium (ISM) emission mechanisms. 

This dust-correlated emission, referred to as anomalous 
microwave emission or microwave excess, represents an addi- 
tional component of the diffuse emission at microwave wave- 
lengths. Most of our early knowledge of this type of emission 
has come from statistical analyses of large areas of sky (de 
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Oliveira-Costa et al. 2002; Lagache 2003; Davies et al. 2006; 
Miville-Deschenes et al. 2008). In recent years, a few detec- 
tions from specific types of Galactic objects have been re- 
ported in the literature, such as Hii regions (Dickinson et al. 
2007, 2009; Todorovic et al. 2010), dark clouds (Finkbeiner 
et al. 2002; Casassus et al. 2006; AMI Consortium: Scaife 
et al. 2009a,b; Dickinson et al. 2010; Vidal et al. 2011) and 
molecular clouds (Watson et al. 2005; Casassus et al. 2008; 
Tibbs et al. 2010). 

Currently, the most popular hypothesis to explain the 
observed microwave excess is that of electric dipole radiation 
from rapidly rotating very small dust grains, also commonly 
referred to as 'spinning dust' (Draine & Lazarian 1998; Ali- 
Haimoud, Hirata & Dickinson 2009). The spinning dust 
model generates a very distinctive peaked spectrum that 
rises between ^ 10-30 GHz, which cannot be easily repro- 
duced by either free- free or synchrotron emission. Alterna- 
tive models have been proposed, including hard synchrotron 
radiation (Bennett et al. 2003) and magnetic dipole radia- 
tion (Draine & Lazarian 1999), however, the spinning dust 
paradigm is the one that best reconciles model predictions 
with observations (e.g. see Planck Collaboration 2011a). 

Although the anomalous emission occurs in the mi- 
crowave regime, the proposed spinning dust hypothesis im- 
plies that the origin of the emission is intrinsically associated 
with interstellar dust. A number of dust models have been 
postulated to explain the infrared (IR) emission in the dif- 
fuse ISM (Desert, Boulanger & Puget 1990; Li & Draine 
2001; Draine k Li 2007; Compiegne et al. 2011). Despite 
differences between models, it is generally accepted that the 
observed IR emission falls into two categories: 1) stochastic 
emission and 2) thermal equilibrium emission. The stochas- 
tic emission is attributed to polycyclic aromatic hydrocar- 
bons (PAHs) as well as to a population of small carbona- 
ceous dust grains, while the thermal equilibrium emission 
is likely to be produced by large carbon and silicate dust 
grains. PAHs, first discovered by Gillett, Forrest & Merrill 
(1973), are large molecules 1 nm in diameter) containing 
up to a few hundred carbon atoms. They are believed to 
be responsible for the broad emission features found in the 
near-IR spectrum (Leger & Puget 1984). The small carbona- 
ceous dust grains are ^2-8 nm in diameter and thought to 
produce the mid-IR continuum emission. Finally, the large 
carbon and silicate grains, from a few nm to a few hun- 
dreds of nm in size, emit continuum radiation at far-IR wave- 
lengths. Current spinning dust models (Ysard & Verstraete 
2010; Hoang et al. 2010; Silsbee et al. 2011) attribute the 
anomalous emission to the smallest dust grains. 

Given the key role played by dust in the understanding 
of the microwave excess, an accurate characterisation of the 
dust properties in regions exhibiting anomalous emission is 
necessary to shed light on the nature of this effect. In this 
work we focus on the Perseus cloud, which to date is the 
most compelling example of anomalous emission associated 
with a specific Galactic region. The Perseus molecular cloud 
complex is a giant molecular cloud chain ^ 30 pc in length, 
containing a total gas mass of ^ 1.3x10^^ located 
at a distance of approximately 250 - 350 pc (Ridge et al. 
2006a). The Perseus cloud is interesting, and hence exten- 
sively studied, not only because of its relative proximity, but 
also because it forms low to intermediate mass stars, thus 
providing a link between the well known low mass star for- 
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Figure 1. Reprocessed Spitzer MIPS 24 /j^rn. map of G159.6-18.5 
overlaid with the VSA 33 GHz contours illustrating the high spa- 
tial correlation that exists between the microwave and IR emis- 
sion. The bright central region of G159.6-18.5 corresponds to the 
BO V star HD 278942, believed to be responsible for the Hii bub- 
ble. The VSA ^ 7 arcmin FWHM synthesized beam is displayed 
in the bottom left-hand corner of the image, and the contours 
correspond to 10, 20, 30, 40, 50, 60, 70, 80 and 90 per cent of the 
peak flux, which is 200 mJy beam"-*^. 



mation in the Taurus molecular cloud and the high mass star 
formation in the Orion molecular cloud. Within the cloud 
is the extended dust structure G 159. 6-1 8. 5. Observations 
performed as part of the COMPLETE Survey (Ridge et al. 
2006a) identified G159.6-18.5 as an expanding Hii bubble 
located on the far side of the cloud. Watson et al. (2005) ob- 
served G159.6-18.5 with the COSMOSOMAS experiment 
on degree angular scales, and found a rising spectrum be- 
tween 11 and 16 GHz. Follow-up observations performed 
by Tibbs et al. (2010) with the Very Small Array (VSA) 
interferometer at 33 GHz on angular scales of ^ 7 - 40 ar- 
cmin, identified five emission features displaying a remark- 
ably high spatial correlation with the mid-IR emission (see 
Fig. 1). The Perseus cloud has also been selected for exten- 
sive investigations of the microwave excess by the Planck 
Collaboration (2011a). By combining the nine Planck fre- 
quency bands in the range 28.5 - 857 GHz with ancillary 
data, and by subtracting the contaminating contributions 
arising from free-free emission, thermal dust emission and 
the CMB, the authors were able to produce an accurate and 
comprehensive spectrum of the anomalous emission on de- 
gree angular scales. 

The layout of the paper is as follows. Section 2 intro- 
duces the Spitzer data used in this analysis and describes 
the data reprocessing. In Section 3 we combine the Spitzer 
data with a dust emission model (dustem) to characterise 
the emission properties of the cloud. In Section 4 we use the 
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results of the dust modelling to characterise the spinning 
dust emission and we present our conclusions in Section 5. 



2 SPITZER DATA 

In order to characterise the dust properties across the 
Perseus cloud and investigate the link with the anomalous 
emission mechanism, we used multi-band IR data from the 
Spitzer Space Telescope, taken as part of the 'From Molecu- 
lar Cores to Planet-Forming Disks' (c2d; Evans et al. 2003) 
Spitzer Legacy Program. The principal aim of the c2d pro- 
gram is to investigate the formation process of stars and 
planets, from prestellar cores to planet-forming disks. To 
achieve this goal, the c2d program involved both photo- 
metric and spectroscopic observations of five nearby molec- 
ular clouds, including the Perseus cloud. The photomet- 
ric data are from the Multi-band Imaging Photometer for 
Spitzer (MIPS; Rieke et al. 2004) and the Infrared Array 
Camera (IRAC; Fazio et al. 2004) while the spectroscopic 
observations were carried out with the Infrared Spectro- 
graph (IRS; Houck et al. 2004). In the present paper, we 
limit our analysis to the MIPS and IRAC data, which pro- 
vide a coverage of the Perseus cloud comparable to the VSA 
33 GHz observations. 

Modelling the dust properties in the diffuse ISM de- 
pends critically on the quality of the IR observations. It 
is therefore crucial that we minimise any systematic effects 
that could potentially bias the analysis. Possible system- 
atics in the Spitzer data include stripes, gaps or artifacts, 
point sources and contamination from zodiacal light. Both 
the point sources and zodiacal light contamination are more 
pronounced at the shorter Spitzer wavelengths. The density 
of sources increases with decreasing wavelength, while the 
zodiacal light peaks around 10 - 20 fim (Fixsen & Dwek 
2002), where its contribution can be comparable to the dif- 
fuse ISM emission.^ For these reasons, we decided to repro- 
cess both the IRAC and MIPS data. Such a reprocessing, in 
particular for the MIPS data, benefited from the techniques 
recently developed by the MIPSGAL (A 24 and 70 /xm Sur- 
vey of the Inner Galactic Disk with MIPS') Spitzer Legacy 
Program team (Carey et al. 2009), which have been opti- 
mised for the treatment of extended emission in highly struc- 
tured background regions such as the Galactic Plane and the 
Perseus cloud. 

We note that, for our reprocessing, we used, respec- 
tively, versions S16 and S14 of the Spitzer Science Cen- 
ter (SSC) pipeline for the MIPS and IRAC Basic Calibrated 
Data (BCDs). The c2d analysis was instead carried out with 
version S13. The reprocessing of the photometric Spitzer 
MIPS and IRAC data is described in Sections 2.1 and 2.2, 
respectively. 



2.1 MIPS Data Reprocessing 

The MIPS instrument observes at three bands centred on 
24, 70 and 160 //m, with angular resolutions of 6, 19 and 
40 arcsec, respectively. Table 1 lists the MIPS data used in 
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Table 1. Summary of the MIPS observations used in the present 
analysis. ^ Units of Right Ascension are hours, minutes and sec- 
onds. ^ Units of Declination are degrees, arcminutes and arcsec- 
onds. 
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AOR Key 




(J2000) 


(J2000) 


Epoch 1 


Epoch 2 


per2 


03:44:30.7 


+32:06:08.1 


5781248 


5787904 


perS 


03:42:34.5 


+31:55:33.0 


5781504 


5788160 


per 4 


03:40:39.0 


+31:37:53.0 


5781760 


5788416 


per5 


03:37:42.0 


+31:16:41.0 


5782016 


5788672 



this analysis. Only the c2d Astronomical Observation Re- 
quests (AORs) that provided coverage of the VSA 33 GHz 
observations were selected. The 24 fim data required mini- 
mal reprocessing, while the 70 and 160 //m data, had to be 
substantially reprocessed. The difference in data quality be- 
tween 24 /im and 70 and 160 /im is caused by the different 
detector materials: the 24 /im detector is a Si As array, while 
the 70 and 160 /xm detectors are Ge:Ga arrays. In particular, 
the responsivity of the Ge:Ga arrays varies significantly with 
the build-up of charge from cosmic ray hits, while the SiAs 
array is less sensitive to this effect (Gordon et al. 2005). 



2.1.1 24 fim Data 

At 24 fim, all the MIPS first frames have a shorter ex- 
posure time. This results in a ^ 10 - 15 per cent sup- 
pression of the signal in the first frame. This effect also 
affects the second and third frames, although to a lesser 
extent, with a 2 per cent suppression. For this reason, 
the first three frames of each observation were discarded. 
This removal had no significant impact on the data as we 
had sufficiently many frames available to compensate. Af- 
ter discarding the first three frames, we subtracted the zo- 
diacal light contribution from each BCD using the model 
by Kelsall et al. (1998), which is implemented in the SSC 
pipeline. This model is derived from observations with the 
Diffuse Infrared Background Experiment (DIRBE) onboard 
the Cosmic Background Explorer (COBE) satellite. Finally, 
an overlap correction routine, developed for the MIPSGAL 
survey (Mizuno et al. 2008), was applied to match the back- 
ground level of all the BCDs. The final mosaic was produced 
by combining the individual BCDs using the SSC MOsaicker 
and Point source Extractor (MOPEX) software package. 



2.1.2 70 iim Data 

For the MIPS 70 /im array, frequent (^ every 2 min) stimula- 
tor fiashes provide a continuous measurement of the respon- 
sivity. However, if a stimulator fiash occurs simultaneously 
on top of a bright region in the sky, it can yield an inaccu- 
rate estimate of the responsivity, and this in turn can have 
adverse effects, such as the production of streaks in the data. 
To remove artifacts due to responsivity drifts in our Ge:Ga 
data, we followed the procedure applied to the MIPSGAL 
70 jim. data as described in Paladini et al. (in preparation). 
The first step of this reprocessing makes use of the Germa- 
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nium Reprocessing Tools (GeRT) software,^ one of the data 
reduction tools provided by the SSC. The GeRT software 
allows the user to run the GeiGa SSC pipeline offline, and 
customise its performance to his/her specific needs. Before 
running the GeRT software on our data, all the stimula- 
tor flash measurements were inspected and any erroneous 
measurements were removed using a moving median filter. 
The GeRT processing was then performed using a piecewise 
linear function to interpolate between the measurements. 
These newly created BCDs were then corrected with a flat- 
fielding routine to remove slowly varying gain fluctuations 
in the detector (see Fig. 2). 

At the same time of the GeRT reprocessing and the flat- 
fielding, it was identified that a thermal anneal took place 
during observations of this region, creating an artificial step 
in the map (see Fig. 2). A thermal anneal occurs approx- 
imately every three hours, raising the temperature of the 
GeiGa detectors to ^ 8 and 4 K for the 70 and 160 /im ar- 
rays, respectively, thereby erasing the accumulated effects of 
cosmic ray hits. To correct for the difference in responsivity 
before and after the anneal, we used IRIS (Improved Repro- 
cessing of the IRAS Survey; Miville-Deschenes & Lagache 
2005) data at 60 /xm. This involved scaling the IRIS 60 /xm 
map to 70 /im, assuming a standard SED for diffuse dust 
emission, and convolving the MIPS 70 /xm map to the 4 ar- 
cmin angular resolution of the IRIS data. After performing 
a colour correction to account for the different bandpasses, 
we compared the IRIS and MIPS maps and found that the 
best match occurred after the thermal anneal had been per- 
formed. By performing five cuts across the step in the map, 
the median offset was found to be 10.3 MJy/sr, and the inac- 
curate BCDs were then corrected for this effect (see Fig. 2). 
As for the MIPS 24 fxm data, the 70 fxm data were then cor- 
rected for zodiacal light, although at this wavelength the ef- 
fect is less significant. An overlap correction was performed, 
and finally, the corrected BCDs, ignoring stimulator flash 
frames, were combined with MOPEX to produce the final 
mosaic. 

2.1.3 160 ^im Data 

The 160 jim. array, like the one at 70 /im, is a Ge:Ga de- 
tector. Hence, the 160 fim. data had to be reprocessed using 
the GeRT software as described in Section 2.1.2. The newly 
produced BCDs did not, however, require flat-flelding, as the 
responsivity drift is in this case less severe than at 70 /im. At 
this wavelength the zodiacal light is negligible (< 2 per cent) 
and therefore its contribution did not have to be subtracted 
from the data. No overlap correction was required. 

The 160 /im detector has a much smaller field of 
view (2.8 arcmin^) compared with the 24 and 70 /im ar- 
rays (29.2 and 13.7 arcmin^, respectively). This, combined 
with the fact that over half of the detector has damaged pix- 
els, results in an incomplete coverage of the map (see Fig. 3). 
To mitigate this effect, we slightly undersampled the image 
by selecting 40 arcsec pixels when producing the final mo- 
saic, allowing MOPEX to interpolate over the gaps in the 
coverage. The map was then regridded from this coarser grid 
onto the native pixel grid (see Fig. 3). 

^ http:/ /ssc. spitzer.caltech.edu/dataanalysistools/tools/gert/ 



Table 2. Summary of the IRAC observations used in the present 
analysis. ^ Units of Right Ascension are hours, minutes and sec- 
onds. ^ Units of Declination are degrees, arcminutes and arcsec- 
onds. 
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2.2 IRAC Data Reprocessing 

The IRAC instrument observes over the wavelength range 
3.0 /im to 10.5 /im, in 4 bands centred on 3.6, 4.5, 5.8 and 
8.0 /im, all with an angular resolution of ~ 2 arcsec. Ta- 
ble 2 lists the IRAC data used in this analysis. Only the c2d 
AORs that provide coverage of the VSA 33 GHz observa- 
tions were selected. IRAC observations are calibrated using 
point sources, but because we are interested in the diffuse 
emission, an extended emission correction was applied to the 
data, following the prescription described in the IRAC Data 
Handbook. All four IRAC bands suffer from well known ar- 
tifacts such as muxbleed, column pulldown and electronic 
banding. To correct for these artifacts, as well as for the 
first frame effect, the IRAC BCDs were processed using con- 
tributed software provided by Sean Carey.^ After performing 
this artifact mitigation, the zodiacal light was subtracted. 
This subtraction was carried out individually on each BCD. 
It should be noted that when the IRAC observations were 
processed by the SSC pipeline, a dark field subtraction was 
performed, in order to correct for dark current and bias off- 
sets. This sky dark observation also contains some contribu- 
tion of zodiacal light, which then must be added back to the 
observations to avoid over-subtraction. To remove large scale 
gradients, a fiat-fielding was performed. Finally an overlap 
correction was applied. As for the MIPS maps, the final mo- 
saic was produced using MOPEX. 

2.3 Point Source Extraction 

All Spitzer bands, with the exception of MIPS 160 /im which 
is characterised by a very low source density, contain a sig- 
nificant number of point sources that require subtracting. 
The source extraction performed at all wavelengths (exclud- 
ing MIPS 160 /im for which no source subtraction was per- 
formed) made use of the Astronomical Point source Extrac- 
tor (APEX) tool, which is part of the SSC's MOPEX soft- 
ware. The APEX source detection and extraction consists of 
median filtering the background and subsequently detecting 
the sources using the Point Response Function (PRE). Fol- 
lowing the recommendation by the SSC, for the MIPS chan- 

^ http:/ /web. ipac.caltech.edu/staff/carey/irac_artifacts/ 
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Figure 2. Spitzer MIPS 70 /im maps: Left: c2d map displaying the vertical striations caused by the variations in responsivity in the 
Ge:Ga detector and the artificial step caused by the thermal anneal; Middle: Our map showing the result of running the GeRT software 
and the fiat-fielding which has removed the striations observed in the c2d map; Right: Our map after correcting for the step caused by 
the thermal anneal. Source subtraction has not been performed on any of these images. 
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Figure 3. Spitzer MIPS 160 /im maps: Left: c2d map displaying the vertical striations caused by the variations in responsivity in the 
Ge:Ga detector and the gaps in the coverage; Middle: Our map showing the result of running the GeRT software which has removed the 
striations observed in the c2d map; Right: Our map after correcting for the lack of coverage using 40 arcsec pixels when creating the 
mosaic with MOPEX. Source subtraction has not been performed on any of these images. 



nels we performed source subtraction on the final mosaics 
(single- frame mode), while for IRAC we first stacked all the 
BCDs and then detected the sources simultaneously in each 
frame (multi-frame mode). In addition, at 24 and 70 fim we 
built a customised PRF by selecting the brightest point-like 
sources in the cloud, while for the IRAC bands we used the 
in-flight IRAC PRFs, recently made available by the SSC.'^ 
For the MIPS 24 //m data, instead of using the stan- 
dard APEX tool, we adopted the MIPSGAL point source 
extraction pipeline, which is optimised for detecting sources 
in regions of highly varying background, such as the Galac- 
tic Plane or the Perseus cloud. The MIPSGAL pipeline is 
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a combination of APEX modules, idl routines and shell 
scripts, and consists of: 1) generating the core PRF by se- 
lecting sources in low background regions and with little or 
no structure; 2) combining this PRF with the wings of the 
theoretical PRF; 3) using the filtered PRF to extract the 
source flux and the unfiltered PRF to subtract the source 
from the input mosaic; 4) applying a scaling factor to correct 
the source flux for the median filtering. More details can be 
found in Shenoy et al. (in preparation). 

2.4 Final Reprocessed Maps 

The IRAC 3.6 and 4.5 fim mosaics are characterised by such 
a high density of sources that, even after performing source 
extraction, it makes them difficult to use to analyse ex- 
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Figure 4. Reprocessed Spitzer IRAC 8 /im image at 40 arcsec 
angular resolution. 



tended emission. In addition to this, the contribution of stel- 
lar emission decreases, while the contribution of PAH emis- 
sion increases, across the four IRAC bands, making IRAC 
band 4 (8 fim) the optimum band to trace the PAH emis- 
sion. Also, the IRAC observations at all four channels were 
performed simultaneously, and due to the positioning of the 
IRAC detectors in the focal plane, the coverage common to 
both IRAC channels 3 and 4 is less than that of the indi- 
vidual channel coverage. For these reasons, we decided to 
ignore IRAC bands 1, 2 and 3 and use only IRAC band 4. 
We note that this channel, together with the three MIPS 
channels, provides coverage of both the stochastic and ther- 
mal equilibrium dust emission. In particular, the stochastic 
emission, produced by PAHs and small carbonaceous dust 
grains, will be dominant at 8 and 24 /xm, while the thermal 
equilibrium dust emission, produced by the large carbon and 
silicate dust grains, will dominate at 160 fxm. The emission 
at 70 /im is a result of both stochastic and thermal equilib- 
rium dust emission (Compiegne et al. 2010). 

The total coverage of our reprocessed maps is 2.3 deg^ 
for the IRAC 8 /xm map, 4.4 deg^ for the MIPS 24 /xm 
map, 3.9 deg^ for the MIPS 70 /xm map and 4.2 deg^ for the 
MIPS 160 /im map. These maps are displayed in Figs. 4, 1, 2 
and 3, respectively. Conservative uncertainties for the flux 
in the Spitzer maps, based on the intrinsic calibration of 
the instruments combined with the uncertainties introduced 
due to the data reprocessing, were estimated to be 5 per 
cent for the 8 and 24 /jui maps, and 15 per cent for the 70 
and 160 /im maps. 



3 MODELLING THE IR DUST EMISSION 

To investigate the interstellar dust we combined the Spitzer 
data discussed in Section 2 with the dust emission model 
DUSTEM (Compiegne et al. 2011). dustem, an updated ver- 
sion of the Desert, Boulanger & Puget (1990) model, al- 
lows the user to parameterise various dust emission and ex- 
tinction properties. The software package incorporates five 
dust grain species: neutral and ionised PAHs (log-normal 
size distribution with centre radius ao = 0.64 nm and width 
a = 0.1), small amorphous carbon grains (ao = 2.0 nm. 



a = 0.35), large amorphous carbon grains (power-law size 
distribution with a minimum radius amin = 4.0 nm, a cut-off 
radius of at = 150 nm, and a spectral index of a = —2.8) and 
large silicate grains {amin = 4.0 nm, at = 200 nm, a = —3.4). 
For further details we refer the reader to Compiegne et 
al. (2011). For this analysis we combined the neutral and 
ionised PAHs and merged the large amorphous carbon and 
large silicate grains, creating three dust grain populations of 
PAHs, Very SmaU Grains (VSGs) and Big Grains (BGs). 

3.1 Constraining the IR Dust Emission 

To constrain the dust properties in the Perseus cloud we 
wrote an SED fitting routine, based on the idl routine MP- 
FIT (Markwardt 2009), that combines the computed dust 
emissivities from dustem with the instrumental bandpass 
filters, and allows us to estimate the flux that would be ob- 
served by an instrument at a given wavelength. Using this 
fitting routine, we would like to constrain the mass abun- 
dances of the PAHs and VSGs relative to the BGs (Ypah and 
Yvsg), the strength of the interstellar radiation field (ISRF) 
relative to the Mathis, Mezger & Panagia (1983) solar neigh- 
bourhood value (xisrf) and the column density of hydro- 
gen (Nh). However, due to a lack of data at wavelengths 
longer than 160 /xm, we do not measure the peak emission of 
the thermal equilibrium dust, which prevents us from recov- 
ering a reasonable estimate of the temperature and thus the 
strength of the ISRF. Consequently, there is a degeneracy 
between the strength of the ISRF and the column density 
of emitting material. 

To circumvent this problem, we make use of the an- 
cillary data from the COMPLETE survey (Ridge et al. 
2006b). This database is complementary to c2d, and in- 
cludes extinction maps and emission maps of atomic and 
molecular line data. In particular, the COMPLETE extinc- 
tion map of Perseus was generated by applying the NICER 
algorithm (Lombardi k Alves 2001) to 2MASS data. Us- 
ing Nh/Av = 1.87x10^^ H cm"^ mag"^ from Bohlin et 
al. (1978), we converted this extinction map into a map 
of the column density of hydrogen, hence breaking the de- 
generacy described above. Since the angular resolution of 
the 2MASS/NICER extinction map is 5 arcmin, this sets a 
lower limit to the angular resolution at which we can con- 
duct our analysis. Taking into account that the VSA obser- 
vations have an angular resolution of 7 arcmin, we degraded 
aU four Spitzer maps, along with the 2MASS/NICER ex- 
tinction map, to this resolution. After regridding them onto 
a 90 arcsec pixel grid, we incorporated the maps within our 
SED fitting routine. 

Having estimated Nh, we keep this parameter fixed in 
the model, and we fit the 8, 24, 70 and 160 /xm data points 
for YpAH, Yvsg and xisrf, resulting in one degree of free- 
dom (dof). We do not constrain the hardness of the ISRF, 
but make the assumption that it is similar to that of the solar 
neighbourhood. We note that if this is not the case, and if the 
ISRF is in fact harder in the Perseus cloud, then this would 
have the effect of decreasing the estimated abundances of 
both the PAHs and VSGs. Nevertheless, we believe that us- 
ing the solar neighbourhood value is a reasonable approach 
in this case, given the lack of high-mass star formation in 
the cloud (Luhman et al. 2003). 

Running dustem and our SED fitting routine on a pixel 
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Figure 5. Maps of the Perseus cloud: Ypah map (Top Left); YvsG map (Top Right); XlSRF map (Middle Left); Nh map (Middle 
Right); T^ust map (Bottom Left); and the VSA 33 GHz map (Bottom Right). The Ypah 5 YysG Xisrf parameter maps are the first 
to be produced for the Perseus cloud. Ypah and YysG cire given relative to the value for the diffuse high galactic latitude medium (see 
Compiegne et al. 2011), while Xisrf is dimensionless. The Nh map was computed from the 2MASS/NICER extinction map and is in 
units of 10^0 H cm . The VSA map is provided to allow a visual comparison, and is overlaid with the thirteen regions listed in Table 3. 
All maps are displayed at 7 arcmin angular resolution. 



by pixel basis resulted in a set of parameter maps (Ypah, 
YvsG and Xisrf) at 7 arcmin angular scales. In addition, for 
each pixel the SED fitting procedure also retrieves a tem- 
perature for each population of grains (PAHs, VSGs and 
BGs). These temperatures are not parameters in the fit- 
ting process, but are computed directly by dustem based 



on both the optical and thermal properties of the grains 
and the results of the fit. More precisely, given an exciting 
ISRF, DUSTEM computes the temperature distribution for 
each size bin of each dust population taking into account 
the multi-photon absorption processes, based on the method 
described by Desert, Boulanger & Shore (1986). Since only 
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the BGs are in thermal equihbrium, we estimate the median 
for the BG size distribution only, and we take this value 
as representative of the dust equilibrium temperature in a 
given pixel (Tdust)- Fig. 5 illustrates the parameter maps 
obtained for Ypah, Yvsg, Xisrf and Tdust- For complete- 
ness, we also show the map of Nh derived from the COM- 
PLETE extinction map. Ypah and Yvsg are given relative 
to the value for the diffuse high galactic latitude medium (see 
Compiegne et al. 2011), Xisrf is dimensionless and Nh is in 
units of 10^° H cm . Remarkably, the maps of Ypah, Yvsg 
and Xisrf are the first of their kind ever generated for the 
Perseus region, providing a comprehensive view of the dust 
temperature, radiation field and abundances of the different 
dust populations across the cloud. 

A careful inspection of these maps reveals that Ypah 
and Yvsg both vary between ^ 2 and 10 throughout 
the bulk of the cloud. Nh ranges between ^ 4x10^^ 
and 2x10^^ H cm~^, and in regions where Nh is 
high 1x10^^ H cm~^), both Ypah and Yvsg appear 
to be decreased, suggesting that grain coagulation may be 
occurring (e.g. Stepnik et al. 2003; Flagey et al. 2009). 
Xisrf ranges from ^ 0.5 to 6, while the Tdust varies be- 
tween ^ 12 to 22 K throughout the cloud. In addition, xisrf 
and Tdust clearly trace the dust shell of G159.6-18.5, with 
both peaking around the open cluster IC 348. Finally, Yvsg 
seems to be decreased in regions where Xisrf is increased, 
which may indicate that the VSGs are being destroyed due 
to the increase in xisrf- 

3.2 Investigating the Robustness of the Dust 
Modelling Analysis 

To cross-check our results and investigate bias effects due 
to the adopted extinction map, we used two independent 
Nh constraints, in addition to the 2MASS/NICER map. 
Firstly we created an Nh map using the gas maps (Hi 
and ^^CO) from the COMPLETE survey. The Hi obser- 
vations have an angular resolution of 7 arcmin and a spec- 
tral resolution of 0.32 km s~^. By integrating the Hi data 
cube, assuming a spin temperature of 135 K (Dickey et al. 
1978), the Hi data were converted into a Hi column density 
map. The ^^CO integrated map has an angular resolution 
of 46 arcsec and it was converted to a H2 column density 
map using the standard relation X = Nh2/W(^^C0), where 
W(^^CO) is the integrated line intensity. We assumed a typ- 
ical X factor value of 2x10^° cm"^ K"^ km"^ s (Strong & 
Mattox et al. 1996). The total Nh map was produced by 
combining both maps (Nhj + 2xNh2) at a common angu- 
lar resolution of 7 arcmin. Secondly, we used an extinction 
map of the region released as part of the c2d data prod- 
ucts at an angular resolution of 5 arcmin. This map was 
generated in a similar manner to the 2MASS/NICER ex- 
tinction map, however, in this case, both near- and mid-IR 
(1.25 - 24 /im) observations were folded in, allowing us to 
probe much higher extinction values (up to Ay ~ 15 mag) 
than our original 2MASS/NICER map. This new extinction 
map was convolved to 7 arcmin. Both these new maps were 
used separately in our SED fitting routine to estimate Nh 
and both gave similar results to those obtained with the 
2MASS/NICER extinction map. Due to the intrinsic vari- 
ations in these three separate Nh constraints, the median 
variation in the parameters across the entire map was com- 
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Figure 6. Scatter plot showing the comparison between the tem- 
perature map obtained in this analysis with dustem and the 
temperature map produced by SGS08. The line of best fit is dis- 
played (dashed line) along with the line y = x (dash-dot line). 



puted to be 5, 6, 13, 23 and 2 per cent for Ypah, Yvsg, 
Xisrf, Nh and Tdust, respectively. The large uncertainty on 
Nh is likely due to variations in the X factor across the 
cloud, as well as to the possible existence of a 'dark gas' 
component - likely associated with a mixture of cold Hi and 
warm H2 - that is not traced by either 21 cm measurements 
or by CO line observations (see, for instance, Planck Collab- 
oration et al. 2011b). Having tested the Nh constraint and 
quantified the uncertainties, we combined them in quadra- 
ture with the fitting errors to estimate the final uncertainties 
on the parameters (see Section 4). 

We also checked the reliability of the dustem analysis 
by comparing our temperature map with temperature maps 
previously obtained for this region. Temperature maps of the 
Perseus cloud have been created by Schnee, Goodman & Sar- 
gent (2008, hereafter SGS08) who used various combinations 
of MIPS 70 and 160 fim and IRIS 60 and 100 fim maps to 
compute colour temperature maps. To compare their tem- 
perature map, and the one produced in this analysis, we 
convolved their 70/160 temperature map to 7 arcmin, re- 
gridded it onto our 90 arcsec pixel grid and trimmed it to 
match the coverage of our temperature map. A scatter plot 
of the two temperature maps is displayed in Fig. 6 and we 
find that the temperature map computed with dustem is 
slightly lower (by approximately 20 per cent) than the tem- 
perature map of SGS08. There may be several reasons to 
explain this discrepancy. As discussed in Section 3.1, the 
temperature computed by dustem is a physical temperature 
based on the optical and thermal properties of the grains, 
while the temperature produced by SGS08 is a colour tem- 
perature based on a grey-body emission model. It is also 
possible that the median temperature we compute for the 
BG size distribution might not be completely representative 
of the effective temperature of the observed spectrum. An- 
other possibility might be due to the fact that SGS08 used 
the original c2d data in their analysis. In this work we have 
completely reprocessed the c2d data, as discussed in Sec- 
tion 2, and this is particularly important for the MIPS 70 
and 160 /im maps due to the variations in the responsiv- 
ity of the Ge:Ga detector. Last but not least, SGS08 might 
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Table 3. Central coordinates of the thirteen 12 arcmin x 12 arcmin regions used in this analysis, and whether or not the region 
corresponds to an anomalous emission feature in the Tibbs et al. (2010) VSA map. Also displayed are the mean parameters computed 
using DUSTEM within each of the thirteen regions. ^ Units of Right Ascension are hours, minutes and seconds. ^ Units of Declination 
are degrees, arcminutes and arcseconds. ^ Abundances relative to the BG abundance. ^ Scaling factor relative to the Mathis, Mezger & 
Panagia (1983) solar neighbourhood radiation field. 



Region 


RA« 


DEG^ 


Anomalous 


YpAH^ 


YvsG^ 


XlSRF^ 


Tdust 


Nh 


xVdof 




(J2000) 


(J2000) 


Emission 




















(K) 


(1020 H cm-2) 




A 


03:44:32 


+32:11:40 


Yes 


2.70 


± 


0.44 


2.14 


± 


0.52 


4.47 


± 


0.88 


21.7 ± 0.7 


89.6 ± 20.6 


0.6 


B 


03:44:56 


+31:53:09 


No 


4.56 


± 


1.46 


4.76 


± 


1.78 


0.68 


± 


0.22 


16.0 ± 9.3 


94.5 ± 21.7 


21.9 


C 


03:44:49 


+31:39:15 


No 


6.01 


± 


2.13 


4.45 


± 


1.89 


0.49 


± 


0.18 


15.1 ± 5.1 


102.0 ± 23.4 


31.6 


D 


03:43:01 


+31:43:27 


Yes 


2.64 


± 


0.47 


2.27 


± 


0.58 


1.44 


± 


0.30 


18.1 ± 1.3 


121.6 ± 28.0 


4.8 


E 


03:43:16 


+32:00:26 


Yes 


2.26 


± 


0.39 


2.23 


± 


0.53 


1.76 


± 


0.36 


18.7 ± 2.5 


123.7 ± 28.4 


5.2 


F 


03:42:07 


+31:54:20 


No 


2.00 


± 


0.34 


2.17 


± 


0.50 


1.42 


± 


0.29 


18.0 ± 1.1 


133.5 ± 30.7 


3.3 


G 


03:41:30 


+31:40:59 


No 


2.23 


± 


0.40 


2.56 


± 


0.61 


1.33 


± 


0.28 


17.8 ± 1.4 


109.6 ± 25.2 


3.4 


H 


03:40:27 


+31:27:34 


Yes 


2.47 


± 


0.42 


2.39 


± 


0.57 


1.81 


± 


0.37 


18.7 ± 1.4 


85.3 ± 19.6 


1.8 


I 


03:39:12 


+32:25:05 


No 


2.72 


± 


0.53 


2.86 


± 


0.74 


1.18 


± 


0.27 


17.5 ± 2.1 


91.8 ± 21.1 


5.1 


J 


03:37:51 


+31:22:07 


Yes 


2.91 


± 


0.56 


2.88 


± 


0.75 


1.51 


± 


0.34 


18.2 ± 1.7 


86.9 ± 20.0 


4.7 


K 


03:37:02 


+31:05:55 


No 


3.84 


± 


0.93 


3.56 


± 


1.11 


0.76 


± 


0.20 


16.2 ± 6.1 


92.1 ± 21.2 


12.5 


L 


03:40:28 


+31:14:05 


Yes 


2.90 


± 


0.49 


2.12 


± 


0.55 


1.92 


± 


0.39 


18.9 ± 1.9 


86.9 ± 20.0 


3.0 


M 


03:39:54 


+31:43:43 


No 


2.80 


± 


0.78 


7.02 


± 


2.21 


1.15 


± 


0.33 


17.4 ± 1.3 


72.8 ± 16.8 


5.5 



have under estimated the contribution of VSGs at 70 /im, 
which would result in a higher dust temperature estimate. 
Nevertheless, the two temperature maps show an overall cor- 
relation which, given that the two maps have been produced 
using significantly different methods, provides robustness to 
our analysis. 



4 CHARACTERISING THE ANOMALOUS 
MICROWAVE EMISSION 

The goal of this paper is to identify any potential correla- 
tions between spatial variations of the dust properties and 
the observed anomalous emission in the Perseus cloud. The 
parameter maps of Ypah, Yvsg, Xisrf, Nh and Tdust pro- 
vide a strong leverage to carry out such an investigation. 
For this purpose, we use the VSA 33 GHz map shown in 
Fig. 5 as a reference. Since the VSA is an interferometer, 
the large scale structures are removed from the map due to 
incomplete sampling of the Fourier plane. However, as dis- 
cussed in detail in Tibbs et al. (2010), one can clearly see 
that, on the spatial scales to which the instrument is indeed 
sentitive, some regions are characterised by the presence of 
anomalous emission, while others are not. We then use this 
information to investigate how the dust parameters in these 
regions compare with the parameters in regions where no 
anomalous emission has been detected. 

For this purpose, thirteen regions were selected in which 
to perform a detailed analysis (see Fig. 5). Each of these 
thirteen regions is 12 arcmin x 12 arcmin, providing cover- 
age throughout the majority of the Perseus cloud. Table 3 
lists the central coordinates of the regions and whether or 
not the region corresponds to an anomalous emission fea- 
ture in the VSA map. Also quoted in Table 3, and plotted 
in Fig. 7, is the mean value of each of the parameters (Ypah, 
Yvsg, Xisrf, Nh and Tdust) within each region. The mean 
value was computed assuming that all the pixels in each re- 



gion are independent, and to assess whether the mean was 
a good statistical indicator in this case, we inspected the 
pixel value distribution for each region, and checked that it 
was approximately gaussian. The uncertainties on the pa- 
rameters were estimated by combining the errors from the 
fit, which were computed from the covariance matrix, with 
the uncertainty due to using the 2MASS/NICER extinction 
map to constrain Nh (see Section 3.2). Additionally, Ta- 
ble 3 lists the mean x^/<^o/ values obtained from the SED 
fitting procedure. The dustem fits appear to be reasonable 
for the majority of the regions, except for regions B, C and 
K {x^/dof> 10). We believe the reason for this is that these 
regions are much colder, therefore the peak of the SED is at 
wavelengths considerably longer than 160 /im. As a conse- 
quence, the current set of data points allow us to only loosely 
constrain the IR emission in these regions. A better repre- 
sentation of their physical dust properties will be achieved 
by incorporating the Herschel SPIRE measurements, ob- 
served as part of the Herschel Gould Belt Survey (Andre et 
al. 2010), which will become publicly available at the end of 
the proprietary period. 

By splitting the thirteen regions according to their lo- 
cation relative to the anomalous emission, a weighted mean 
was computed for each of the parameters listed in Table 3. 
The results of this analysis are illustrated in Table 4. Both 
these results and Fig. 7 suggest a significant increase in 
Xisrf, and hence Tdust, in regions of anomalous emission. In 
addition, we find that there is a slight decrease in the abun- 
dance of VSGs in regions of anomalous emission, however, 
this is only a 1.2cr result and is therefore not statistically sig- 
nificant. We also find that there is no significant variation 
in the abundances of the PAHs or the column density of hy- 
drogen between regions with and without anomalous emis- 
sion. Interestingly, despite an enhancement in Xisrf (and 
Tdust) in regions of anomalous emission, there does not ap- 
pear to be any obvious correlation between the strength of 
the anomalous emission observed with the VSA (represented 
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Figure 7. Plots displaying the computed parameters listed in Table 3 for all thirteen regions, along with the mean value {dashed line). 
Regions with {filled squares) and without {open squares) anomalous emission are displayed, with the size of the squares representing the 
strength of the anomalous emission observed with the VSA (Tibbs et al. 2010). These plots show that there appears to be a clear increase 
in both xisRF and T^ust in regions of anomalous emission, but that there does not appear to be any obvious correlation between the 
strength of the anomalous emission and the parameters computed in this analysis. 



by the size of the squares in Fig. 7) and these parameters. 
This is also true for the other parameters. This, and the 
fact that some of the results in Table 4 are only tentative, 
implies that anomalous emission is a very complex process 
to parameterise and is clearly dependent on a combination 
of multiple parameters, some of which we were not able to 



constrain in this analysis (e.g. ionisation state and size dis- 
tribution of the PAHs). 

The recent investigation by the Planck Collaboration 
(2011a) finds that the anomalous emission in the Perseus 
cloud can be accurately modelled by spinning dust emission 
arising from two gas components: a high density molecular 
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Table 4. Results of performing a weighted mean analysis on 
the parameters in Table 3. These results show that for regions 
with anomalous emission there is a significant enhancement in 
the strength of the ISRF and the dust temperature. 



Parameter 



Anomalous 
Emission 



No Anomalous 
Emission 



Significance 



YpAH 

YvsG 

XISRF 

Nh 

Tdust 



2.60 ± 0.18 
2.29 ± 0.23 
1.74 ± 0.15 
94.3 dz 9.0 
20.2 di 0.5 



2.44 ± 0.21 
2.75 ± 0.31 
0.87 d= 0.09 
93.8 ± 8.3 
17.7 lb 0.7 



0.6(7 
1.2cr 
4.9o- 
O.Oo- 
3.1o- 



gas component and a low density atomic gas component. 
The results of our analysis can also be explained in terms of 
the spinning dust paradigm. Since we find that there is no 
apparent significant variation in the abundances of PAHs 
and VSGs between regions with and without anomalous 
emission, the grain excitation mechanism will be the most 
important factor in producing spinning dust emission. The 
enhancement in Xisrf in regions of anomalous emission will 
result in a higher density of photons which, in turn, will 
increase the photon - grain interactions, hence producing 
more excitation that could cause the dust grains to spin. 

With the VSA, Tibbs et al. (2010) found that the mea- 
surable microwave - IR correlation is restricted to the dense 
knots located along the shell of the Hii bubble G 159.6- 
18.5 and the open cluster IC 348. It is known that the ra- 
diation pressure within expanding Hii regions sweeps up 
the surrounding dust forming a Photo Dissociation Re- 
gion (PDR; Hollenbach & Tielens 1999). This indicates that 
the knots of anomalous emission observed with the VSA 
along the shell are originating from the PDR surrounding 
the Hii gas, excited by the central BO V star HD 278942. 
As an open cluster, we expect an increase in the strength 
of the ISRF observed towards IC 348 due to the presence 
of young stellar objects that are known to reside there (e.g. 
J0rgensen et al. 2008), thus possibly providing the source 
of excitation required to produce the observed anomalous 
emission in IC 348. Both of these scenarios are consistent 
with our results of finding an enhancement in xisrf in re- 
gions with anomalous emission, indicating that xisrf, and 
hence photon - grain interactions, play a vital role in the 
excitation of the dust grains. 

The Planck Collaboration (2011a) results show in par- 
ticular that the bulk of the anomalous emission is associated 
with the dense molecular gas component, and in such dense 
environments, gas - grain collisions are the predominant ex- 
citation mechanism of the dust grains. However, it must be 
remembered that the VSA is an interferometer, and as such 
it automatically filters out large scale structures, implying 
that the anomalous emission detected by the VSA is differ- 
ent to that detected by Planck. In fact, Tibbs et al. (2010) 
calculated that the VSA only observes approximately 10 per 
cent of the total 33 GHz emission while Planck is sensitive 
to the entire cloud. Given that the star formation efficiency 
in the Perseus cloud is only 2.6 per cent (J0rgensen et al. 
2008), the bulk of the cloud would appear to be quiescent - a 
result which is consistent with the anomalous emission de- 



tected by Planck being due to a dense molecular gas com- 
ponent. Moreover, the difference in the angular responses of 
the VSA and Planck would seem to explain why xisrf ap- 
pears to play an important role for the anomalous emission 
observed with the VSA, but not for the Planck observations 
at degree angular scales. 

Anomalous emission has been observed in a number 
of Hii regions (e.g. Dickinson et al. 2007; Todorovic et al. 
2010), however, the observations lack the resolution to re- 
solve their structure, thereby making it difficult to determine 
the exact location from which the emission is originating. 
G159.6-18.5 just happens to be large enough that the VSA 
observations can resolve the Hii region interior and identify 
that the anomalous emission is indeed produced in the PDR. 
The idea of anomalous emission originating from a PDR is 
not new, as - for instance - anomalous emission has previ- 
ously been observed in a PDR in the p Ophiuchi molecular 
cloud (Casassus et al. 2008). 

To further constrain the spinning dust emission model, 
we decided to combine the parameter values obtained with 
DUSTEM with SPDUST (Ali-Haimoud, Hirata & Dickinson 
2009; Silsbee et al. 2011), a software package that com- 
putes a spinning dust spectrum for a set of user defined 
input parameters. Here we focus on xisrf as we have al- 
ready shown that this parameter is enhanced in regions 
of anomalous emission. Taking the value of Xisrf for each 
of the thirteen regions listed in Table 3 and adopting all 
other SPDUST parameters used by the Planck Collaboration 
(2011a) for the Perseus cloud (uh = 250 cm"^, Tgas = 40 K, 
Uh+Zuh = 1.12x10"'^, nc+/nH = 1x10""^, 2nH2/nH = 1), 
we estimated the corresponding spinning dust curves. Fur- 
thermore, we used our values of Xisrf in combination with 
the idealised environmental parameters for other phases of 
the ISM derived by Draine & Lazarian (1998), to compute 
spinning dust curves for a Dark Cloud (DC), Cold Neutral 
Medium (CNM), Warm Neutral Medium (WNM), Warm 
Ionised Medium (WIM), Reflection Nebula (RN) and PDR. 
In this work we used the updated version of spdust as 
described by Silsbee et al. (2011), which assumes that the 
grains are randomly orientated with respect to their angular 
momentum axis. The computed spinning dust curves for one 
of the regions with anomalous emission (region E) are dis- 
played in Fig. 8 along with the spinning dust emissivity of 
the anomalous emission observed with the VSA (Tibbs et al. 
2010). Region E was chosen as it corresponds to the location 
of the peak anomalous emission detected by the VSA. 

From Fig. 8 it is clear that the spinning dust models for 
each of the phases of the ISM are not consistent with the 
VSA measurement, with the model for the WNM represent- 
ing the best agreement. Although we are using the values 
of Xisrf that we have derived in this analysis, we believe 
that the reason for this lack of agreement is because the 
other idealised values for the environments are not physi- 
cally representative of the region. For instance, the idealised 
value of Xisrf for the Draine & Lazarian (1998) PDR model 
is 3x10^, which is clearly not consistent with the values of 
Xisrf obtained in this work (see Table 3), suggesting that 
the other idealised environmental values are also not appro- 
priate. The result for the other regions of anomalous emis- 
sion is similar, implying that to fit the observations we need 
more constraints on the other parameters within spdust, 
rather than using idealised values. 



12 C. T. Tihhs et al 




10 100 1000 

V (GHz) 



Figure 8. Spinning dust curves computed using spdust with the 
XlSRF values constrained with dustem for region E and ideahsed 
values for all other parameters. Also plotted is the spinning dust 
emissivity computed from the VSA observations. This shows that 
the modelled spinning dust curves are not consistent with the 
VSA observations, and that xisrf alone provides very little con- 
straint on the computed spinning dust emission. This highlights 
the need for more constraint on the other spdust parameters. 
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Figure 9. Spinning dust curves computed with spdust for all 
thirteen regions using the xisrf values constrained with dustem 
and the idealised parameters for the WNM (Top). The differ- 
ence in these spinning dust curves relative to region A (Bottom). 
This shows clearly, that even in regions with no anomalous emis- 
sion, the modelled spinning dust emissivity is similar to regions in 
which we detect anomalous emission. This implies that these mod- 
elled curves cannot explain the observed emissivity variations, and 
once again highlights the complexity of this emission mechanism. 



We also find that even in regions where we do not detect 
any anomalous emission, the modelled spinning dust emis- 
sivity is similar to regions where we do detect anomalous 
emission. This idea is shown in Fig. 9 where the modelled 
spinning dust emissivity for the WNM is plotted for all thir- 
teen regions: remarkably, the computed spinning dust curves 
have very little dependence on Xisrf- This highlights that, 



although Xisrf was found to play a significant role in the 
excitation process of the dust grains, the strength of the 
ISRF alone is not enough to accurately constrain the spin- 
ning dust emission and distinguish between regions with and 
without anomalous emission. This, in turn, once again em- 
phasises the complex nature of this emission mechanism and 
demonstrates that other parameters, not yet constrained by 
this analysis, may play a bigger role. 



5 CONCLUSIONS 

This work represents the first attempt to characterise the 
anomalous emission by using IR observations to constrain 
the dust properties. By reprocessing the original c2d Spitzer 
photometric observations of the Perseus cloud we have pro- 
duced maps of a substantially improved quality. Using these 
reprocessed maps and the 2MASS/NICER extinction map 
from the COMPLETE survey, in conjunction with the dust 
emission model dustem, we have been able to constrain the 
abundances of the PAHs and VSGs relative to the BGs, 
the strength of the ISRF relative to the Mathis, Mezger & 
Panagia (1983) solar neighbourhood value and the thermal 
equilibrium temperature of the BGs, resulting in parameter 
maps of the region at 7 arcmin angular scales. The Ypah, 
YvsG and Xisrf parameter maps represent the first of their 
kind for the Perseus cloud. 

We found that there is a significant increase in xisrf, 
and consequently Tdust, in regions of anomalous emission. 
This result can be interpreted in terms of the spinning dust 
hypothesis, with the anomalous emission possibly being pro- 
duced by electric dipole emission from spinning dust grains 
excited by an increase in photon - grain interactions. This 
result also suggests that the spinning dust excitation mecha- 
nism plays a vital role in producing the anomalous emission 
observed with the VSA and as such the emission observed 
with the VSA is likely originating from the PDR surround- 
ing the Hii gas in G159.6-18.5, and the open cluster IC 348. 

The results of our analysis also suggest that, due to 
the difference in the observed angular scales, the anomalous 
emission detected with the VSA is different to that detected 
by the Planck Collaboration (2011a), which would appear 
to explain why we find that Xisrf plays an important role 
for the anomalous emission observed with the VSA, but not 
for the Planck observations at degree angular scales. 

In this work we demonstrate the complex nature of the 
anomalous emission. With more IR observations becoming 
available (e.g. Herschel data), it will be possible to constrain 
additional dust properties to even higher angular resolution. 
Combining this with future high resolution microwave obser- 
vations will help to identify, in even more detail, the physi- 
cal conditions in which anomalous emission originates. This 
type of analysis creates a new perspective in the field of 
anomalous emission studies, and represents a powerful new 
tool for constraining spinning dust models. 
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